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FOREWORD 


The work summarized in this report was performed 
under contract NAS8-27980, Space Shuttle Engine 
Program, Coca 4B Reactivation Review. 


ABSTRACT 

The investigation and analysis of a high-pressure 
hydrogen facility compressor is chronicled, and a 
life prediction based on fracture mechanics is 
presented. Crack growth rates in SA 105 Gr II 
steel are developed for the condition of sustained 
loading using a hypothesis of hydrogen embrittle- 
ment associated with plastic zone reverse yielding. 
The resultant formula is compared with test data 
obtained from laboratory specimens. 
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SUMMARY CONCLUSIONS 


The following conclusions have been reached as a result of the subject analysis 

1. The safe operating life at 15,000-psi hydrogen maximum operating 
pressure is in excess of 100 cycles. 

2. The safe operating life at 10,000-psi hydrogen maximum operating 
pressure is in excess of 460 cycles. 

3. The equivalent damage in terms of life cycles caused by partial 
excursions in pressure has been estimated. 

4. Conditions hazardous to personnel exist at hydrogen pressures above 
5000 psi during the first half of the safe operating life (before 
the first 230 cycles at 10,000 psi). 

5. Conditions hazardous to personnel exist at all hydrogen pressures 
after the first half of the safe operating life (after the first 
230 cycles at 10,000 psi). 
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BACKGROUND 


In February of 1974, this office was requested to assess the feasibility of per- 
forming a fracture mechanics-based life prediction analysis of the high-pressure 
hydrogen facility compressor at Coca 4B test site. It was concluded that such an 
analysis could be performed to indicate qualitatively the criticality of the com- 
pressor for the intended service. Crack growth rate data (cyclic) were available 
for a material similar to the one from which the compressor was constructed, and 
the compressor structural details were largely known. The 4-inch-thick, welded, 
23-inch-diameter, high-pressure spheres were selected to be the most critical 
components, and life prediction analyses were performed. In recognition that 
sustained- load flaw growth was not accounted for in the fracture data (cyclic), 
a life factor of 2 was used to reduce the nominal prediction. It was concluded 
that the compressor life was indeed critical and significantly less than the 
required service life. 

The next several months were spent largely in a technical interchange with the 
customer. Most of the effort was spent in performing additional calculations 
and supplying background information, both at customer request. Also during this 
period, additional crack growth data were obtained for cycle periods of 1 and 10 
seconds/cycle. The fracture specimens were cut from forgings of the same material 
as the subject pressure vessels, and the cyclic data confirmed the values used in 
previous analyses. At the end of this period, there was clearly a technical dis- 
agreement between this office and the customer regarding the criticality of the 
compressor. It was decided to obtain an independent assessment by Dr. P. Paris 
of Del Research Corporation, an internationally recognized consultant. 

Dr, Paris was on-site at Canoga on 1 June. This morning was spent in briefings 
intended to provide Dr, Paris with all of the basic information available con- 
cerning the compressor system, the materials, and the hydrogen-environment frac- 
ture data. Dr. H. Marcus of the Rockwell International Science Center also con- 
tributed in the discussions. This office intentionally remained out of the 
briefings to ensure the integrity of Dr. Paris* independent assessment. The mo§t 
significant event of the morning briefings was the presentation of sustained- load 
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crack growth data (not previously shown) by Dr. W, Chandler.. When Dr. Paris gave 
his preliminary assessment in the afternoon, he focused on the sustained- load 
data and suggested that failure might occur in as few as 10 cycles. In his docu- 
mented report of 7 June, Dr. Paris expressed essentially the same opinion. 

This office did not accept the above-mentioned sustained- load data at face value, 
and pursued an inquiry to determine the conditions of testing. It was found that 
all loadings which resulted in a significant amount of crack growth had been pre- 
ceded by at least 1000 cycles of fatigue cycling (in hydrogen) at the same stress 
% 

intensity level. This observation led to a new theory which seemed to explain 
the existing data and suggested a means of quantifying sustained- load crack growth 
for service conditions (where no fatigue cycling occurs between sustained loads). 

Meetings were held on 12 and 13 June at the Brown University (Providence, Rhode 
Island) office of Dr. Paris to review the sustained- load data and to discuss other 
aspects of the life prediction analyses. At the conclusion of these meetings. 

Dr. Paris was largely in agreement with this office, and he documented his new 
position with an addendum to his previous report. By mutual agreement, it was 
recommended that additional growth rate data be obtained at loading periods of 
100 and 1000 seconds/ cycle. 
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BASIC ASSUMPTIONS 

The following assumptions form the basis for fracture mechanics life prediction of 
the hydrogen compressor system: 

1. The most critical system component is the 23- inch-diameter, 4-inch- 
thick, welded, high-pressure (15,000-psi] hydrogen sphere. 

2. The most critical and likely location of an undetected crack for a 
welded sphere is in the weld bead, on the inner surface (exposed to 
hydrogen), and oriented normal to the weld pass direction. 

3. The crack shape aspect ratio is 1/2, i.ei, the surface length (major 
diameter of a semi-ellipse) is twice the crack depth. 

4. The high fracture toughness of the subject material precludes the use 
of proof test logic for determining the maximum flaw size after proof 
test. 

5. The maximum probable flaw depth after proof test (initial flaw size in 
service) is 1.00 inch. 

6. The weld material stress, normal to the crack, during proof loading is 
0.80 times the material ultimate tensile strength. This stress repre- 
sents a combined effect of welding residual stresses, applied proof 
stresses, and weld material yielding. 

7. The weld material stress, normal to the crack, during operation is 
linearly proportional to the ratio of the operating pressure divided 
by the proof pressure. 

8. The *’Irwin" part-through crack equation adequately predicts the applied 
stress intensity for crack depths up to and including one-half the wall 
thickness. Stress intensity values for deeper cracks are unpredictable. 

9. The sustained-load crack growth due to hydrogen- embrittlement effects 
is one-tenth of the plastic zone size calculated from the maximum load 
reached during the cycle. 
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10. The mechanical fatigue component of the total growth per cycle is small 
relative to the hydrogen- embrittlement effect. 

11. The above assumptions are sufficiently conservative to preclude the 
necessity for additional safety factors in life prediction calculations. 
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CRACK GROWTH RATES 


During the course of the subject investigation considerable effort was required to 
establish the crack growth rates that would occur during service operating con- 
ditions. Real-time sustained-load holding periods of 6 hours or more precluded 
the gathering of very much laboratory data under service conditions. However, 
an analysis of available data obtained at much smaller loading times (5 minutes/ 
cycle and 1 hour/ cycle) suggested a hydrogen- embrittlement fracture mechanism 
that could be proven Coi* disproven) with relative ease in the laboratory, using 
load holding times of 100 seconds/cycle and 1000 seconds/cycle. It was recognized 
from smooth bar fatigue test results that hydrogen- embrittlement effects were 
directly related to plastic flow in the material. Extending this knowledge to 
material conditions near the crack front, the following fracture mechanism was 
postulated: 

1. Hydrogen diffusion is largely determined by plastic strain. 

2. Sustained- load growth occurs where hydrogen is diffused. 

3. Crack growth arrest occurs where hydrogen is not diffused. 

4. After many cycles in fatigue, hydrogen is throughout plastic zone. 

a. Sustained load grows crack through plastic zone. 

b. Crack arrest occurs at end of plastic zone. 

5. Without fatigue cycles, hydrogen diffuses only a small depth. 

a. Depth assumed same as material seeing full hysteresis loop. 

b. Crack growth per cycle equals hydrogen depth. 

6. Hydrogen diffusion is also time-dependent . 

Dr. E. F. Cain recognized that if hydrogen was not present during plastic material 
strain, then embrittlement should not occur. Tests using a helium environment 
during loading and unloading (during plastic straining), but also subjecting the 
specimen to high-pressure hydrogen during the load-holding period (no plastic 
strain), resulted in no measurable crack growth. Although not conclusive, the 
helium- hydrogen testing tended to confirm the postulated fracture mechanism. 
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An attempt to quantify item 5a above was based on theoretical estimates of reversed 
yielding (during unloading) near the crack front. Estimates ranged from one-tenth 
to one-quarter of the tension plastic zone size. It was reasoned that material 
which did not undergo reversed yielding would have much less plastic strain through 
the cyclic hysteresis loop, much less hydrogen diffused through the material, and 
much less crack growth under sustained load. It will be shown in the following 
pages that the assumption of a one-tenth plastic zone size growth per cycle is 
conservative for this material in a 15,000-psi hydrogen environment. 

At this time, a complete report of material testing is not available for reference 
so a summary of the pertinent data is presented in Table 1. Two stress intensity 
ranges, 20 and 40 ksi ^inch, were selected on the basis of previously obtained 
cyclic data and specimen limitations. The primary variable was the sustained- load 
hold time, with the end objective of determining crack growth during a 6- hour 
duration. It can be seen in the test results that hold time did not have a signifi- 
cant effect at the lower stress intensity. A trend did exist at 40 ksi >/inch, but 
it was not sufficiently definitive to be used for crack growth prediction. As 
expected, the time- dependent growth rate (yin. /sec) was significantly reduced as 
the cycle period was extended. The incremental change between the 10-sec/cycle 
growth and the 100-sec/ cycle growth was 1.22 yin./sec, whereas the time rate of 
growth between 100 sec/cycle and 1000 sec/cycle was 0.27 yin./sec (taking the 
data at face value) . These nianbers support the hypothesis that the growth/ cycle 
approaches an asymptote as the cycle period becomes very large. However, the 
meager amount of data precludes calculation of the asymptote value. The assumed 
crack growth based on plastic zone size is more than three times any of the measured 
values, so it is believed that the time-dependent factor has been adequately 
accounted for. 

Also shown in Table 1 are two data values for ramp loading times of 100 seconds. 

It was noted that the hydrogen compressor system required periods on this order 
to reach maximum pressure, whereas the other data were acquired using a ramp rate 
of 0.5 second. On the theory that plastic strain (and therefore hydrogen diffusion) 
occurs during loading and unloading, it was expected that the ramp rate. could greatly 
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TABLE 1. CEACK GROWTH RATE DATA SUMMARY 
SA 105 STEEL IN 15,000-PSI HYDROGEN 


AK, 

ksi y inch 

R 

Ramp 

Time, 

seconds 

Hold 

Time, 

seconds 

DADN, 

pin. /cycle 

20 

0. 

1 

0. 

5 

9 

120 







9 

220 







99 

81 







999 

200 

20 

0. 

1 

100 

100 

132 

40 

0. 

1 

0. 

5 

9 

440 







9 

480 







99 

570 







999 

810 

40 

0, 

1 

100 

100 

790 


influence the resulting sustained-load growth rate. The data show that, at a 
40-ksi y / inch stress-intensity range, the effect is significant but not over- 
whelming, These data are included in the comparison with the assumed growth 
rate discussed below. 

Figure 1 shows a plot of the data from Table 1 along with the assumed crack growth 
rate of one-tenth the plastic zone size (based on 58,0 ksi material yield strength). 
The dashed line that seems to fit the test data was obtained by dividing the assumed 
growth rate by a factor of 4, hence the similarity in curve shapes. The uppermost 
data point at 20 ksi ^Tich is about one-third of the asstimed rate, and at 40 ksi 
v/inch the ratio is even better. It is assumed therefore that the assumed growth 
rate is adequately conservative. 

In all of the above discussions the crack growth rates have been related to the 
parameter of stress intensity, always on the basis of 15,000-psi hydrogen pressure. 
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Figure 1. Comparison of Assumed Growth Rate With Test Data 
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Obviously, in a pressure vessel the stress intensity must be a reflection of the 
hydrogen pressure. Therefore, a given crack size will have a high hydrogen pres- 
sure associated with a high stress intensity and a low hydrogen pressure associated 
with a low stress intensity. The practical limits of this investigation prevent 
a full study of the pressure effect on growth rate so the analysis must be based 
on the worst condition, i.e., the 15,000-psi hydrogen environment. It is believed 
that these growth rates introduce an additional conservatism in the analytical 
life predictions for pressures less than 15,000 psi. 
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CYCLE LIFE FOR 10,000-PSI PRESSURE 


After reviewing the structural characteristics of the various components of the 
hydrogen compressor system, it seemed evident that the second-stage discharge 
pulse quieter spheres were significantly more critical than any of the other com- 
ponents. Therefore, the compressor system life is determined by the life of these 
spheres. The two major factors that determine the spheres to be critical are: 

(1) the 4-inch- thick welds that are very likely to contain crack-like defects, 
and (2) the virtual impossibility of adequately inspecting the welded joint to 
detect defects. The relatively low strength and high fracture toughness of the 
material could allow an improbably large crack to survive proof test without de- 
tection, so a judgment probability estimate was used rather than proof test logic 
in determining the maximum possible flaw size. Also, consideration was given to 
the size of the weld heat-affected zone in estimating the crack shape. The frac- 
ture mechanics calculations below are based on an initial flaw that is semi-circu- 
lar in shape, oriented normal to the direction of the weld pass, 1.00 inch deep, 
and exposed to the vessel's internal hydrogen pressure. The basic information 
listed below provides a brief structural description of the second-stage discharge 
pulse quieter spheres: 

Sphere Diameter: 23.1 inches OD 

Wall Thickness: 4,3 inches 

Construction: Forged Hemispheres Welded 

Material: ASME SA 105 Gr II Steel 

= 58.0 ksi, = 79,4 ksi 

Proof Pressure: 26,250 psi 

Design Operating Pressure: 15,000 psi 

The Irwin stress intensity equation for part-through cracks is used as a basis 
for fracture mechanics calculations. This is a relatively crude approximation for 
stress intensity because it does not account for the compound curvature of the 
pressure vessel walls but rather assumes the crack to be in a flat plate. Also, 
the equation needs an additional correction factor when the crack depth exceeds 
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half the plate thickness, but because of the compound curvature this will not be 
attempted. Hopefully, the conservatism of stopping crack growth analysis at half 
the wall thickness (2.1-inch crack depth) will offset the unconservatisms inherent 
in applying a flat plate equation to a compound- curvature problem. 

The weld nugget solidifies at a time when large thermal gradients exist across 
the weld joint. As the joint cools after welding, the thermal gradients vanish 
and are reflected as strain gradients, thereby causing high residual tensile stress 
in the weld region. The process is, in fact, a casting process and as in castings 
it is not unusual to create shrink cracks when the stress/ strain level becomes too 
high. Assuming that this type of defect exists in the spheres, then it follows 
that the residual stress (in the direction of welding) is on the order of the 
material yield strength. Adding the proof pressurization stress to the residual 
stress must cause a significant amount of plastic strain, with the stress level 
approaching the material ultimate strength. For this analysis it is assumed that 
the weld region stress at proof pressure is 80 percent of the material ultimate 
strength . 


a . = 0.80 X 79.4 = 63.5 ksi 
proof 

For all pressurizations in service the stresses will be lower and linearly 
proportional . 

*^operating ^proof ^^operating^^proof^ 

For 10,000-psi operating pressure: 


Oio QQQ = 63.5 (10,000/26,250) == 24.2 ksi 
Conservatively calculating the plastic zone size on the basis of a high ratio of 
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where 


K = stress intensity 


= material yield strength 


K = 1.12 ct/WQ = 1.32 a/T 

(a/2c = 0.5, o/0„, = 1.0) 

K. = 1.32 X 24.2 v/^ = 31.9v/a 


operating 


y operating 




51.9y/a r 
58.0 / 


= 0.0161a 


DADN = crack growth rate = ” 0.00161a 


Table 2 uses the above result to arrive at the number of cycles required to grow 
an initial crack 1.00 inch deep to the half wall thickness depth of 2.10 inches, 
assuming a constant crack growth rate over each 0.10-inch increment. Rounding 
the resultant value to two significant figures arrives at a predicted life of 460 
cycles for the critical component of the hydrogen compressor system. 

TABLE 2. SECOND-STAGE DISCHARGE PULSE QUIETER 
CYCLE LIFE FOR 10,000-PSI MAXIMUM PRESSURE 


a 

avg 

DADN = ^ 

A XT ^21 

^ " DADN 

1.05 

0.00169 

59.2 

1.15 

0.00185 

54.0 

1.25 

0.00201 

49.7 

1.35 

0.00217 

46.0 

1.45 

0.00233 

42.9 

1.55 

0.00249 

40.1 

1.65 

0.00265 

37.7 

1.75 

0.00282 

35.5 

1.85 

0.00298 

33.6 

1.95 

0.00314 

31.9 

2.05 

0.00330 

30.3 

Jan='461 
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ACCUMULATED DAMAGE CALCULATION 
FOR 

HIGH-PRESSURE HYDROGEN SERVICE 


The procedure described below is intended to be used for calculating the damage 
accumulation caused by service operation of the high-pressure hydrogen compres- 
sor at SSFL. The actual service operation of the compressor will not always be 
to the same maximum pressure. Also, the pressure excursion will not always be 
a simple variation from zero to maximum to zero. A cyclic damage equivalence 
calculation is necessary to have an accounting procedure for limiting the pres- 
surization in service. The 10,000-psi hydrogen maximum operating pressure has 
been used as a reference basis, and all damage calculations are expressed as 
equivalent 10,000-psi cycles. The operating life is limited to 460 equivalent 
10,000-psi cycles. 

Step 1, Divide the pres sure- time history of the compressor system into a series 
of block loadings. 

Each pressurization, from zero to one or more pressures to zero, constitutes one 
block of loading. The damage accumulation caused by each block of loading is 
determined by the pressure excursions that occur within the block. Figure 2 
illustrates the division of an example pressure-time history into loading blocks 
Note that the number of loading blocks is established by the number of times 
that the pressure returns to zero, and is independent of the higher pressure 
excursions. 

Step 2, Divide each block loading into a set of pressure excursion cycles which 
may be superimposed to reconstruct the block loading. Disregard holding times 
and pressure variations of 5 percent or less as pressure excursions. Pressure 
excursions from to Pj^ to P are equivalent to pressure excursions from P^ to 

Figure 3(a) illustrates equalities to be used in defining pressure excursions. 
Note that when the hold times and the order of pressure levels are disregarded, 
the pressure excursions of 3(a) are equal for the purpose of damage calculations 
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Figure 3(b) applies the results of Figure 3 (a) and the principle of superposi- 
tion to divide the given loading block into two pressure excursion cycles. 
Figure 3(c) further illustrates the same principles. 

Referring back to Fig, 2, it can be seen that block 1 consists of one pressure 
excursion cycle from 0 to Pj to 0. Block 2 consists of one pressure excursion 
cycle from 0 to to 0. Block 3 consists of two pressure excursion cycles: 

0 to to 0, and to ?2 to P^* 

Step 2, CaZoutate the equvoaZent oycti.a damage caused by each pressure cycle 
excursion using the following formula: 

N equivalent oyelee = x ^ 

where: ksij maximum hydrogen pressure during the excursion 

^min ^ minimum hydrogen pressure during the excursion 

Figure 4 provides example pressure cycle excursions that can be used to demon- 
strate the procedure of Step 3, The N equivalent cycles of excursion cycle A 
are calculated below: 

P = 5, P . =0 

max * mm 

N = 5 X (5 - 0)^ X 10'^ = 0.125 

The equivalent cyclic damage of all the pressure excursion cycles of Fig. 3 are 
listed in the summary comprising Table 3. 

Step 4, Sum the equivalent cyclic damage caused by all of the loading blocks 
and the pressure excursion cycles within the loading blocks. 

Figure 5 provides an example of three loading blocks to illustrate the proce- 
dure of Step 4. The pressure excursion cycles were chosen to correspond with 
those of Fig. 4. Table 4 shows the summary calculations leading to the total 
number of equivalent damage cycles caused by the three loading blocks. 
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PRESSURE 



Figure 2. Pressure Loading Blocks 
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Figure 3. Pressure Excursion Cycles 


PRESSURE, KSI 



Figure 4. Example Pressure Cycle Excursions 
for Damage Calculation 
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PRESSURE, KSI 





TABLE 3, EQUIVALENT CYCLIC DAMAGE OF EXAMPLES 
SHOWN IN FIG. 4 


Excursion Cycle 

P 

max 

P . 

min 

N 

A 

5 

0 

0.125 

B 

10 

0 

1.000 

C 

15 

0 

3.375 

D 

15 

10 

0.375 

E 

15 

5 

1.500 

F 

10 

5 

0.250 


TABLE 4. EQUIVALENT CYCLIC DAMAGE OF 
LOADING BLOCKS SHOWN IN FIG, 5 


Block 

Cycle 

P 

max 

P . 
min 

N 

1 

A 

5 

0 

0.125 

2 


B 

10 

0 




F 

10 

5 

0.250 

■ 

3 


C 

15 

0 

3.375 



D 

15 

10 

0.375 



E 

15 

5 

1.500 

Total 

Equivalent Damage Cycles = 

6.625 
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DESIGN LIFE 


The design life requirement for the total program is estimated to be as many as 
2000 operating cycles. An exact number is not required because there is little 
doubt that a material change is required as soon as possible for the critical com- 
pressor parts. The new material is apparently not embrittled by hydrogen, and 
the service life will not be limited after the change is made. However, testing 
requirements demand that a number of pressurizations be performed during the 
interim. 

The design life requirement for the part of the program prior to compressor rework 
(materials change) has not been established. Verbal information indicates that 
adequate life is available to support planned testing. 
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HAZARD TO PERSONNEL 


Failure of a high-pressure hydrogen sphere will most likely result in high-velocity 
steel fragments permeating a large surrounding region. Obviously, great care must 
be taken to have no personnel present when this might occur. Fracture mechanics 
technology has provided valuable information regarding the hazardous conditions 
leading to failure. 

When a compressor system is never operated above 5000-psi hydrogen pressure, there 
is little likelihood that an initial crack has grown significantly during service 
operation. Therefore, there is a relatively small hazard when personnel are present 

When a compressor system has been operated above 500G-psi hydrogen pressure 
(10,000 psi or 15,000 psi), there is a strong possibility that significant damage 
(increase in crack size) was caused by the higher pressures. This negates the 
rationale for concluding that operation at 5000 psi or less is always safe. Also, 
it should be noted that failure on a succeeding cycle may occur at a much lower 
level than previously sustained. This has been proven by both fracture analysis 
and testing. The degree of hazard for 5000-psi operation is determined by how 
much damage was caused by the higher pressures; the more high-pressure cycles, the 
greater degree of hazard. It is recommended that the compressor not be operated 
(at any pressure) with personnel present after half of the fracture mechanics pre- 
dicted life has occurred in service. 
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APPENDIX A 


COMPUTER PROGRAM FOR CALCULATING EQUIVALENT DAMAGE 

A computer program has been devised which automatically converts input data of 
real-time pressure history into equivalent numbers of 10,000-psi cycles, adding 
the incremental damage to the number of accumulated cycles at the data start point 
The computer program thereby provides the total accumulated damage in terms that 
can be compared to equivalent 10,000-psi cycle life limitations. 

In addition to applying the criteria that defines "loading blocks" and "pressure 
excursion cycles," the computer program utilizes the technique of range pairing 
to reconstruct each loading block in terms of equivalent pressure excursion cycles 
A listing of the computer program is shown in Fig. A-l(a) and A-l(b). 

The following information constitutes instructions for use on a General Electric 
Computer terminal. 

As input, the H-440 FORTRAN program requests the following: 

1. Compressor ID (from 1 to 6 characters) 

2. Accumulated cyclic damage (prior) 

3. Run ID (from 1 to 6 characters) 

4. Data (enter pressures (psi) for run) 

Repeat Steps 3 and 4 for a maximum of 500 pressures. To terminate input, press 
carrier retuim key when program requests run ID. 

The program produces the following output: 

1. Equivalent cyclic damage for each cycle 

2. Incremental cyclic damage (211.) 

3. Prior cyclic damage (Input 2.) 

4. Accumulated cyclic damage (22, and 3.) 

Figure A- 2 shows an example pressure history used to test the various program 
functions. Figure A-3 shows the output obtained when the example was input. 

RSS-8597 

23 





\ 

I 


OlOOSNDM 

lOOOC OAMAGli CALCULATION ■* HIGH PRESSUHE HYDROGEN SERVICE 

iOlOC FOR G* A* VrtOMAN BY J* R* WHITNEY AUGUST 1974 

1020C 

1030 DIMENSION D( SOO > # P< 50 ) 

1040 DATA ND/500/* NP/5O/*0/" '*/ 

lObUC 

1060 5 PRINT* "ENTER COMPRESSOR ID" 

lUIO READ 10>C 

I0»0 10 F0RMAT(A6) 

1090 PRINT lb*C 

I09S lb FORMAT!" ENTER ACCUMULATED CYCLIC DAMAGE FOR COMPRESSOR ”*A6) 

1100 READ#TNE 

1110 SN£=0 

1120 N=1 

1130 20 NNsN 

U40 PRINT* "RUN ID?" 

1 1 SO READ 10*D(1*N) 

1160 IF(D!1*N) .EU. B) go TO 300 

1 1 70 DO 30 I = I * NP 

1180 30 P<1>=«999. 

1190 PRINT* "DATA?" 

1200 REAO*P 

1210 DO 40 MM=l*NP 

1220 1F(P<MM) .£0. -999.) GO TO SO 

1230 40 CONTINUE 

1240 MM = NP-»I 

1250 50 MM«MM- 1 

1260 IF<P(MM) .EO. 0.) GO TO SO 

1265 MM»MM+1 

1266 P<MM)=sO 

1270 Ml=2 

1 280 BN> 1 . 

1290 100 DO 110 «M«M1*MM 

1300 1F(P<M) -EO. 0.) GO TO 120 

1310 110 CONTINUE 

1320 M=MM+l 

1330 120 IF<M .N£* MM > M=M- 1 

1340 Il=Ml-l 

1350 1=11 

1360 200 1=1+1 

1370 210 1F<I .GE» M) GO TO 230 

1380 T=(P!l+l >-P< I) )/<P<l 1-P<1-1 > > 

1390 1F(T .LI* 0.) GO TO 200 

1400 DO 220 J=1*MM 

1410 220 P( J)=P( J+1 > 

1420 MM=MM*l 

1430 M=M-| 

1440 GO TO 210 

1450 230 K=0 

1455 IF<M .£0* MM> M*M- 1 

1460 DO 240 1=M1*M 


Figure A-l(a). Computer Listing of Accumulated Damage Calculation Program 
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1470 1FCP<1) *GE. GO TO 240 

U«0 T=P<I) 

MVO P<U=PU-1) 

ISOO P(1-1)=T 

IblO K=l 

lb20 240 COi>iTINUE 

1S30 IF<K .£0. 1) G0 T0 230 

lb40 CN-0. 

1545 M2s(ll+M)/2 

1546 IF<C .£Q« B) PRINTI 1 # I I »M# <P< I)# 1 = 11 >M > 

1547 II F0RrtAr<2I4/<8F8.0>) 

1550 1MAX=M 

1 560 DO 250 l=I 1#M2 

1570 lFtP<l) -EQ* PlIMAX)) GO T0 250 

1575 CN=CN+1 • 

1580 D< WN+l >=D( I#N) 

15y0 0<2»N)=8N 

1600 0(3»N>=CN 

1610 D(4,N)=PC1MAX) 

1620 D<5»N)=P<1> 

1630 D<6*N)=PC IMAX>4<P< 1MAX)-P< I ) )**2* 1 4E- 1 2 

1640 SM£=SNE4‘DC6^N) 

1645 N=N+1 

1650 2b0 1MAX=IMAX-1 
1660 BN=BN>|4 

1670 IF<M»1 »EQ4 MM) G0 T0 20 

1680 Ml=M+2 

1690 G0 T0 too 

1700C 

1710 300 PRINT 310*C 

1720 310 FORMAT!/////** DAMAGE CALCULATION F0R COMPRESSOR **#46// 
1730& •• RON 8LK CYC PMAX PMIN NE*') 

1735 N=N-l 

1740 00 320 J=1#N 

1745 PRINT 330* <OCl*J)# 1=1*6) 

1750 IF<D<l*J> .NE. D<1*J+I)> PRINT* 

1755 320 CONTINUE 

1760 330 FORMAT! 1X*A6*2F5.0*2F8.0*F9 .3) 

1770 ANC=TN£+SNE 

1780 PRINT 340* SNE* TNE* ANE 

1790 340 FORMAT! /8X* "INCREMENTAL CYCLIC DAMAGE"* F9. 3 
1800* //14X*"PRI0R cyclic 0AMAGE"*F9.3 

1810* //8X» "ACCUMULATED CYCLIC DAMAGE"* F9 . 3) 

1820 PRINT 350 

1830 350 FORMAT! //////) 

1835 IF!C .EO. B) GO TO 5 

1840 STOP 

1850 END 


Figure A-l(b). Computer Listing of Accumulated Damage Calculation Program 
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10 


RUN NO. 1 


RUN NO. 2 


RUN NO. 3 


RUN NO. 4 


Figure A- 2 
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10 10 



. Pressure History Used for Example Accumulated 
Damage Calculation 
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RUN 


DAMAG£ 10J40 NR T/S AUG 20# 197A 


£NT£R C0MPR£SSOR lO 
7 TEST 

ENTER ACCUMULATED CTCLIC DAMAGE E0R C0MPRESS0R TEST 
7 0 

RUN 10? 

7 1 


data? 

7 0 SQOO SOOO 10000 

NUN 107 
7 2 

DATA 7 

7 0 10000 bOOO bOOO 

RUN 107 
7 3 

data? 

7 0 10000 bOOO bOOO 10000 

RUN 107 
7 4 


I A r 

7 0 bOOO 2b00 7bOO 5000 10000 5000 75000 RbOO*-*-*-*-*-*- 2500 3000 

RUN 107 
7 


Figure A-3Ca). Example of Accumulated Damage Calculation Computer Output 
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OAMAce CALCULATiaN F0R C0MPHESS0R TEST 


RUN 

BLK 

CYC 

PMAX 

PM IN 

NE 

1 

1 * 

1 * 

10000. 

0. 

1 .000 

2 

1 • 

1 • 

10000. 

0. 

1 .000 

3 

1 * 

1 • 

10000. 

0. 

1 .000 

3 

1 « 

2* 

10000. 

5000. 

0.250 

4 

1 . 

1 • 

10000. 

0. 

1 .000 

4 

1 . 

£2# 

7S00. 

2500. 

0.187 

4 

1 . 

3* 

7500. 

2500. 

0.187 


INCREMENTAL 

CYCLIC 

DAMAGE 

4*625 



PRIOR 

CYCLIC 

DAMAGE 

0.000 


ACCUMULATED 

CYCLIC 

damage 

4.625 


ST0P 

RUNNING TIMES 3.1 SECS 1/0 TIME t 7.9 SECS 
0K 


BYE 

0FF AT 10144 

Figure A-3(b). Example of Accumulated Damage Calculation Computer Output 
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